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Abstract 

The bright sohton in optical fiber is generally investigated via its spatial evo- 
lution in the time domain, where its waveform is considered in many studies. 
To be consistent with the well-established picture of the dynamics of solitons 
in other systems, in this letter, we propose it is helpful to study the temporal 
evolution of the bright soliton by examining its waveshape propagating along 
the space coordinate axis. We develop a singular theory. Equations governing 
the evolution of the parameters of the bright soliton in the slow time and the 
radiated field are explicitly formulated for the first time. In addition, localized 
modes are found to appear. 

PACS (numbers): 03.40.Kf, 52.35.Mw, 42.65.Tg 

Owing to the promising application to long distance soliton-based communication and 
the great fundamental interest of physics of the process involved, solitary waves and solitons 
in the nonlinear monomode optical fiber have received intensive studies in recent years . 
The generalized propagation equation of optical field in the fiber takes the form 
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iu'^, + ikiu[, — -k2u[,^, + a \u'\ u' = ieP' [u] (1) 

in which x' represents the propagation distance, t' the time and u' the complex field envelope. 
Usually, eP' [u], including linear loss, high-order dispersion and other nonlinear effects, is 
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assumed to be small and treated as perturbations to place emphasis on important phenomena 
of the bright and dark solitons in the fiber 0. In the region of anomalous group- velocity 
dispersion (GVD), by introducing the retarded time T' = f — kix' = t' — x' /vg, Eq. (1) is 
normalized as 

iu'^„ + + \u"f u" = leP" [u] (2) 

in terms of T = T'/Tq, x" = x'/Ld = x' \k2\ /Tq and u" = a/ | A;2 1 / aT^u' Customarily, 
Eq. (2) is referred to as optical nonlinear Schrodinger equation (NLSE), and its unperturbed 
version supports distortionless propagation of a type of solitary wave called the bright or 
temporal soliton 0. 

Generally, waves travelling along the x-axis at speed v are expressible as functions of 
{x—vt). A wave t) may be thought of as formed from the shape f\Q by the substitution 
C = (x — ft), or else as built from the time signal hij) by the substitution r = (t — x/v). 
Here, the function f{() with f{x) = F{x,0) characterizes the "waveshape", and h{T) with 
h{t) = F{0,t) depicts the "waveform" Resulting pictures from the two standpoints 
for the wave F{x,t) are that the "waveshape" changes and propagates along the x-axis 
as time elapses and the "waveform" distorts versus the retarded time r as the distance x 
keep increasing. These actually presents two different point of views for the visualization of 
scenario of soliton under perturbations. 

The bright soliton propagating in the fiber governed by Eq. (2) was typically investigated 
by interchanging the roles of the retarded time T and the space x" and defining an "initial- 
value" problem, or equivalently by directly treating the space x" as the evolution coordinate 
and defining a boundary-value problem. Accordingly, the aspect of waveform of the bright 
soliton was taken into consideration and studies could benefit from the direct application 
of the celebrated frameworks developed by Zakharov and Shabat (ZS) and Ablowitz-Kaup- 
Newell-Segur (AKNS). Nevertheless, to avoid complication of the ZS and AKNS schemes, 
other elaborate approaches were developed in the framework of direct expansion as well . 

In contrast to the studies of bright soliton, the aspect of waveshape is extensively ex- 
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amined in other soliton problems with perturbations 0, including envelope soliton of the 
integrable cubic NLSE in water and other applications QJl^] . Although results for the under- 



standing of the waveform of bright soliton have been achieved, a natural question, how the 
waveshape evolves in the real time or what the dynamics of bright soliton is, is inevitable to 
arise. To answer this question, the corresponding mathematical model is essentially different 
from the one investigated in previous theories, and is also intractable in the ZS and AKNS 
schemes. Consequently, a new theoretical challenge turns up. In this letter, we introduce 
our theory for the subject. 

Let's start from the dimensionless form of Eq. (1) in the anomalous dispersion regime 
of the fiber 

iux + iut + -Utt + \uf' u = ieP [u] (3) 

where t = t'/to, x = x'/l = x' \k2\ /t^, u = a/ | A;2 1 / at^u' and to = 1^21 /ki = \dvg/duj\ . 
Obviously, instead of the usual Tq that is determined by the width of the input waveform 
in existing theories, a characteristic time, namely to that is determined by the working 
wavelength and nature of the fiber, is used in the normalization. Formally, Eq. (3) differs 
from Eq. (2) or the normal form of optical NLSE only by an additional term iut due to 
invalidation of the retarded time, but the essential difference lies in that the time t here 
must be treated as the evolution coordinate and an initial-value problem is consequently 
defined, since the waveshape of bright soliton is to be taken into account. In the absence of 
perturbations, the bright soliton admitted by Eq. (3) is given by 

Usoi{x,t) = 2?7sech2?7(2C + - ^2(\ l) ^ ~ ^'^ 

X exp{-t[2ie -V' + C)x - 2Ct - e,]} (4) 

provided the initial waveshape is in the form 

Usoiix, 0) = 2r/sech277(2C + l)[x - x'] 

xexp{-t[2{C-r]^ + 0x-ei]}. (5) 



From the general point of view of solitons under perturbations, if the perturbations turn 
on, the bright sohton with a starting state of Eq. (5) to propagate in the fiber governed 
by Eq. (3) can not be described by Eq. (4), it undergoes a slow change in the time via 
the variation of its parameters. Moreover, other wave modes come on to appear To 
characterize the picture, we introduce a slow time scale ti = et and assume that the solution 
of Eq. (3) is of the form 

u{t, z, ti) = [2r//(^) + ev{t, z, t^)\e-''^^^'-^ (6) 

where /(z) = sech^;, z = 2ri{2(+ l){x — e^'^x ^ x') ^ = {Kz — e^^Oq — 6i). Also, we need 
further to assume that r], x, x'? ^o? are dependent directly on ti. Obviously, z is 
the coordinate variable in the reference frame tied up to the bright soliton. Here, if we take 
t, z and ti in place of t and x as new independent variables, the derivatives with respect to 
time and space in Eq. (3) are thus replaced by 

-e2,(2C + l)xi| + .J^ (7) 

and 

g-^-M2C+l)g; (8) 

where Xh = A is defined. Introducing Eqs. (6)- (8) into Eq. (3), we transform from 
the laboratory frame into the soliton's one and get two equations from 0(1) and 0{e), 
respectively. Examining the zeroth-order equation for 0(1), we derive A = {2( + 1)~^, 
6oti = = 2{('^ + rf)K and K = {C"^ — ff + QAt]^^. By virtue of these relations, we simplify 
the first-order equation for 0{e) as 

^Vtt + iA~^Vt - 2r]v^t + '^rfv^^ + Srfh'^{z)v + Arfh'^{z)v* - 2rfv = R{z) (9) 

where asterisk * denotes the complex conjugate. And the "source term" |]TT[ R{z) = Rr{z) + 
iRi{z) is given by 



and 



Rr = - lm{Pe'') - Ar]{27]t, + vHh)Mz) 

- 27^[A-\r]-\t, - 2Avt,) - 2{2r}ACt, + Vt,)]M^) 
+ 16r]^A~\[^<pl{z) - 8r/[2r/ACi, + Vtr]z<Pl{z) (10a) 

- 2[A~\, + Ar,\r\t, - Ar/iJ]<^i(z) 

- Arf[e,,, + 2^A-'^,K + A-^x;j^2(^) (10b) 

where 0i(z) = sechz, 02(^) = 2;sechz, f^i^z) = sech2(l — ztanhz), (p2{z) = sechztanh^ 
are defined for simplicity and later use. Expectably, a fresh equation comes out after the 
linearization. Here, we should note that although the basic idea of the present linearization is 



a natural extension of the normal scheme of multiple scale expansion the implementation 
of the idea in handling such soliton problems as of the second order derivative with respect 
to time is original. As usual, extra freedoms for the purpose of preventing the occurrence of 
secular terms are introduced and included in the source term. Taking advantage of Laplace 
transform to solve Eq. (9) yields 

^s'^v + iA'^sv -2'qsvz + 2rfv^^ + Srf'h^{z)v + Arfh^{z)v* - 2rf'v = s~^R{z) (11) 

where v stands for the Laplace transform of v. Putting v = vi + if 2 and v = vi + iv2 = 
{wi + iw2)e^ , we derive from the real and imaginary parts of Eq. (11) 

swi + 2r]'^ALiW2 = s~^ARie~^ (12a) 
SW2 - 27fAL2Wi = -s'^ARre'^ (12b) 

where two Hermitian operators Li =d^/dz^ + (2sech^2; — 1) and L2 =d^/d2;^ + (6sech^z — 1) 
are defined. 

To solve Eq. (12) by virtue of eigen-expansion, a complete set of basis is needed. Consid- 
ering the homogeneous counterpart of Eq. (12), we derive the following eigen- value problem 



Licj) = Xip (13a) 
L2ip = X<p. (13b) 

Now, if we define a non-Hermitian operator H = L2L1, then the corresponding adjoint 
operator is = L1L2. Using the operator L2 to act on both sides of Eq. (13a) and then 
the Li on Eq. (13b) gives 

H(f) = X^cj) (14a) 
= AV- (14b) 

Eigenstates of operators H and H"^ are composed of a continuous spectrum with eigen- 
value A = — (/c^ + 1) and doubly degenerated discrete states with eigenvalue A = 0, re- 
spectively. Under the definition of inner product in the Hilbert space, Their eigenstates 
(f) = {(f){z,k),(f)i{z),(j)2{z)} and ip = {ip{z,k),{pi{z),{p2{z)} turn out to be a biorthogonal 
basis (BB) with the completeness relation 

r" + 00 

(f){z, k)^*{z\ k)dk + Mz)Mz') + Mz)Mz') = S{z - z') (15) 



where 



(j)(z,k) = ^= (l-2iktanhz- k'^)e'^' (16) 



and 
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^{z, k) = ^ ^ (1 - 2sech2z - 2ik tanh z - k^Y^^ (17) 

' ^(F + 1)^ ' ^ ' 



represent the continuous spectrum and 0i(^), 02 (^)5 V^i(^)) V^2(^) that are defined above 
stand for the discrete states. BB is popular in the studies of non-Hermitian Hamiltonian 
problems |T^. With the set of BB, we can expand the solutions of Eq. (12) as 



"+00 

U7i(t,2;, ti) = / t^i (t, A;,ti)v9(z, A;)dA;+ Wii (t,ti)v9i(2;)+ ti?i2 (t,ti)v52(-2) (18a) 
-00 
"+00 

W2{t,z,tx)= I W2 it,k,ti)(p{z,k)dk+ W21 it,ti)(pi{z)+ W22 {t,ti)(j)2iz). (18b) 

-00 



Introducing Eq. (18) into Eq. (12) and solving by means of orthogonality of the basis, we 
derive Wi and W2. Thus, v 1 and V2 are determined from the inverse Laplace transformation. 
Some terms directly proportional to t and t"^ are found to appear in vi and V2, they are 
non-physical and called secular terms. But if we require 

/+00 
Ri{z)M^)d^-Q (19) 
-00 

/+00 r+oo 
Ri{z)(j)2{z)dz + 2riA / Rr{z)zip2{z)dz ^ (20) 
-00 J —00 

/■+OC 
Rr{z)(p2{z)dz = (21) 
-00 

/+00 f+00 
Rr{z)(pi{z)dz + 2riA Ri{z)z(j)i{z)dz = 0, (22) 
-00 J —00 

those terms vanish and we then get the final solution 

J —^(sin (3) Ri{z')(f)*{z', k)ip{z, k)dz'dk 

/-j-oo r-j-oo 1 
/ -—(l-cosp)RJz')^*(z',k)^(z,k)dz'dk 

- ^R,{z')z'M^'W^i{z) 

/+00 * r+oo a2 

-Ri{z')z'M^')dz'+ / —Rr{z')z"ip2{z')dz']ip2{z) 
-00 277 2 



-00 ^'1 d —00 

f+00 A /• + CO 



and 



/H-oo A r+co 
—R,{z')(t>2{z')dz'+ / K^Rr{z')z'^2{z)dz']z^2{z) (23) 
-00 J-oo 

/H-OO p-\-oo I 
/ {smP)Rr{z')<p*{z\ k)(t){z, k)dz'dk 

/-^rOO r-\-oo -1 
/ -cos/?}i?i(^')</'*(-2',^)</'(-2,^)d-2'dA; 
-00 J-oo ^^7 

+ [/ -Rr{z')z'^,{z')dz' + / — i?,(y)^'Vl(^')d^>l(^) 
J-oo J-oo ^ 

/+00 A ;'+oo 
R^{z')^,{z')dz' + / A2i?,(^')^>i(^')d^>'/'i(^) 
■00 J-oo 



-00 

f + OO 



+ / —Rr{z)z'^2{z')dz'U^) (24) 



where we define /3 = 2rfK\{t — '•^^^^•* ). It is noteworthy that localized modes turn out 
to appear in the solution, which is essentially different from the envelope soliton of the 



integrable cubic NLSE with the first-order temporal derivative From our viewpoint, 

the localized modes here is a kind of internal modes whose occurrence is acknowledged to 
be intrinsic for nonintegrable models, for instance, the 0^ model |T5[|. Although it is not 
sufficient to conclude that the model we consider here is nonintegrable, the corresponding 
Lax representation is really difficult to find. Returning to the restriction condition imposed 
on the solution, we indicate that they can be satisfied by the extra freedoms we introduce 
in advance. In fact, they result in a sequence of novel equations 



A 



+00 



r]t^ = ^-^ I Re(Pe^^)sech;zd;z (25) 
^ -'-00 



A f^°° 

Ct, = / Im(Pe*^)tanhzsechzd;z (26) 

2 J ~OD 



V[A-2 + lr/2]x;^ 



/+00 r+oo 
Re{Pe''^)zsechzdz - 2r]A / lm{Pe''^)ztaBhzsechzdz (27) 
-00 J —00 

/+00 r+co 
Im(Pe^^)sechz(l - ztanh2)d^ - 277A / Re(Pe^^)2sechzdz, (28) 
-00 J —00 

which govern the dynamic evolution of bright soliton in the time. In accordance with the 

usual definition of the width w = l/2rj{2Q + 6), we can derive a useful equation 



r" + 00 

,2„,,3 



= 4:7] w / Im(Pe* ) tanh 2;sech2;d2; 
J —00 

/+00 
Re(Pe^^)sechzd2. (29) 
-00 



Now, we generate some specific results by examining two cases. At first, we consider 
the linear loss given by P [u] = —aiu. This perturbation leads to Re(Pe*^) = — 2?7Q;isech2;. 
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From Eq. (25), we compute rjt^ = —2air]A, and then we obtain i] = r^oe"^"^^*^ = 7706"^*^"^^* 
by integration. In this case, A remains constant, thus, the propagation distance of a fixed 
point of the sohton is calculated by a; = At. As a result, we can write 

which recovers a well-known result in previous theories 

Secondly, we give a brief study of the perturbation P [u] = —ia2ud\uf/dt accounting 
for the Raman effect. Using Eq. (7), we derive Im(Pe*^) = — 3277^02 tanh2;sech^2, which has 
influence on the soliton's width and velocity. By Eq. (29), we get Wt-^ = —8a2{2i])^w^/15, 
integrating this equation yields 

w = wo[l + —a2{2r]fwlti]-^ , 

which exhibits that the soliton is narrowed under this effect. As well, we can derive that 
the velocity decreases, obeying 

A = Ao[l + ^a2(2#A2ti]-i 

Under the picture of waveform, the width is depicted differently, and the velocity of dynamic 
sense can not be deflned. 

In conclusion, we think that the waveshape presents a more transparent picture of directly 
physical significance than the waveform, especially in the study of soliton under perturba- 
tions. Hence, we believe that our theory is a nontrivial and necessary alternative for the 
subject. Moreover, the mathematical development in this paper is distinct and normal, its 
idea is helpful for the study of other soliton problems as well. 
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